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ABSTRACT 



The interstellar medium (ISM) in galaxies is directly affected by the mass 
^ and energy outffows originating in regions of star formation. Magnetic fields are 

O cin essential ingredient of the ISM, but their connection to the gaseous medium 

and its evolution remains poorly understood. Here we present the detection of a 
gradient in Faraday rotation measure (RM), co- located with a hole in the neutral 
I hydrogen (H l) distribution in the disk of the nearby spiral galaxy NGC 6946. 

The gas kinematics in the same location show evidence for infall of cold gas. The 
^ combined characteristics of this feature point to a substantial vertical displace- 

'— ' ment of the initially plane-parallel ordered magnetic field, driven by a localized 

star formation event. This reveals how the large-scale magnetic field pattern in 
galaxy disks is directly influenced by internal energetic phenomena. Conversely, 
^ magnetic fields are observed to be an important ingredient in disk-halo interac- 

^ tions, as predicted in MHD simulations. Turbulent magnetic fields at smaller 

\Q spatial scales than the observed RM gradient will also be carried from the disk 

^ and provide a mechanism for the dynamo process to amplify the ordered mag- 

I netic field without quenching. We discuss the observational biases, and suggest 

>■ that this is a common feature of star forming galaxies with active disk-halo flows. 

X 

5^ Subject headings: Galaxies: evolution — Galaxies: individuahNGC 6946 — 

Galaxies: ISM — Galaxies: magnetic fields — Galaxies: structure 



1. Introduction 



The disk-halo interface is crucial for understanding the evolution of spiral galaxies. Star 
formation drives vertical hot gas flows in a process which is described by the galactic fountain 
( Shapiro fc Field||1976 Bregman||1980 ) or chimney model ( Norman fc Ikeuchi||1989 ). In both 
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models, the outflowing hot gas eventually condenses into cool clouds that fall back down 
onto the disk on timescales of tens of Myr. These returning clouds are thought to make 
up at least a part of the population of high and intermediate velocity clouds (HVCs and 
IVCs). The returning cool gas flow must be interacting with the ambient disk-halo medium, 



as is indicated by the HVC and IVC morphology (Briins et al. 2000 Heitsch & Putman 



2009 ) as well as the kinematics of thick disks in external galaxies and the Milky Way ( Heald 



et al. 2007 Fraternali & Binney 2008 Kalberla & Kerp 2009). Magnetic flelds appear to be 
essential in reproducing the structure and evolution of the multiphase ISM (e.g., de Avillez 



& Breitschwerdt 2005 Hill et al. 2012) and the longevity of infalling clouds (Santillan et al. 



& Breitschwerdt 


2005 


2004 


)• 



NGC 6946 is a nearby {D = 6.8 Mpc; Karachentsev et al. 2000) grand-design spiral 



galaxy that has been studied at a broad range of wavelengths. It has a substantial integrated 
star formation rate (2.8M0yr~^; Calzetti et al. 2010), together with a large number of 



catalogued holes in the Hi density distribution (Boomsma et al. 2008) that have sizes, 
masses, and energetics consistent with formation by localized pockets of star formation 
activity. For many of these features, the stars which provided the required energy are 
not distinctly visible, but this is typical of H I holes (e.g.. Brinks &; Bajaja 1986 Bagetakos 



et al. 



2011 



Warren et al. 



2011). The properties of the magnetic fleld (B) in NGC 6946 
have also been studied in detail, showing that its energy density is comparable to other ISM 
components, and perhaps even dominates at large galactocentric radii (Beck 2007). The 
large-scale magnetic fleld has been probed using polarized synchrotron radiation and its 
Faraday rotation measure (RM oc J n^B ■ d/ where is the thermal electron density and / is 
the line of sight), both of which trace the ordered fleld. Faraday RM is a sensitive tracer of 
the line of sight magnetic fleld, and is robustly measured using the RM Synthesis technique 



(Brentjens & de Bruyn 2005; Heald et al. 2009). As in other spirals, NGC 6946 has an 



axisymmetric spiral pattern in the disk combined with a quadrupolar poloidal component 



that becomes dominant at large vertical distances from the disk (Braun et al. 2010). High- 



quality Hi line data (Boomsma et al. 2008) and radio continuum polarization data from 



a combined dataset spanning 1300 — 1432 and 1631 — 1763 MHz (Heald et al. 2009) are 
available in the literature. 



2. A magnetized Hi bubble 



The image of RM (Heald et al. 2009) corresponding to the widespread diffuse polarized 



synchrotron radiation in the disk of NGC 6946 is displayed in Figure [T} and reveals a re- 
markable feature at (a, (5)j2ooo.o = (20'*35"^18^2, 60°06'18"). There, a clear RM gradient is 
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co-located with one of the H I holes (see also Figure |2]) previously identified in the literature 
(Boomsma et al. 2008). The RM at that location ranges from 18.8 — 57.0radm~^. The 



midpoint of the RM gradient is close to the typical RM value in the surrounding area of 
the disk, about 40radm~^. The RM gradient is on a scale of 17", or a linear size of 0.6 
kpc (somewhat smaller than the Hi hole itself, which is 0.9 kpc in diameter from rim to 
rim). The position angle of the gradient is about 225°, similar to the local orientation of 
the ordered magnetic field as can be clearly seen in Figure [TJ By measuring the typical 
fiuctuations in the RM values on the same angular scale throughout the rest of the galaxy, 
this particular feature is found to be significant (i.e., distinguishable from a noise feature) 
at the 4(7 level. Here, a includes both noise and contributions from real fluctuations in B 
and rif. on 0.6 kpc scales throughout the star forming disk. Very few RM fluctuations with 
similar magnitude are found, and the others are preferentially at the edge of the polarized 
disk with low confldence. We can also assess the probability of a chance association of this 
resolved RM feature with an H I hole. This is done by computing the density of H I hole cen- 
ters within the regions with well-detected diffuse polarized emission, taking into account the 
quoted positional uncertainties of 4" x 4" ( Boomsnia||2007 ) . In this way, we estimate that the 



chance of a spurious spatial association of the RM feature and the H I hole is approximately 
0.5%. 

The RM gradient is centered at the typical RM value in the surrounding disk, making 
it unlikely to be caused by fluctuations in rig. The RM gradient is, on the other hand, 
naturally explained by a signiflcant vertical deviation of the local ordered magnetic fleld 
direction - the fleld orientation pointing toward the observer on one side, and away on the 
other. The observed RM pattern originates on the front side of the disk because turbulence 



in the midplane depolarizes radiation from the backside (e.g., Braun et al. 2010). Since RM 
is deflned to be a positive quantity for a LOS magnetic fleld pointing toward the observer, 
the observed pattern implies that the fleld is oriented toward the observer on the side marked 
'A' in Figure [T| and away on side 'B'. Thus the ordered magnetic fleld is directed inward 
(along the spiral pattern toward the center of the galaxy). This agrees with the conclusion 



drawn from the much larger-scale RM gradient already observed across the entire disk ( Beck 
2007|. 



3. Discussion 



The coincidence of the RM gradient with the H I hole strongly suggests that this mag- 
netic structure is caused by the energetic process that led to the formation of the H I hole 
itself. The H I properties of the hole indicate that it is a fairly young (r ~ 20 Myr; Boomsma 
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Fig. 1. — RM image of NGC 6946 from Heald et al. (2009). Magnetic field vectors are shown 
with hnes. The Hi hole catalogued by Boomsma et al. (2008) and found to be colocated 
with the RM gradient is indicated with a black ellipse. 
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Fig. 2. — Images of NGC 6946. All panels are presented on the same angular scale and 



indicate the locations of the Hi holes catalogued by Boomsma et al. (2008) with black 



ellipses. Insets show the immediate vicinity of the feature described in the text. Top left: 
H I column density. The blue line shows the slice used to create the PV diagram in Figure [3j 
Top right: RM image. Bottom left: optical image (created from DSS-II B,R,IR plates). 



Bottom-right: GALEX Nearby Galaxy Survey (NGS; |Gil de Paz etaL][2007| FUV-NUV 
false color image. 



2007) but well-defined structure (albeit with low Hi column density contrast; see Figure [2]). 
A position-velocity (PV) diagram along a slice through the hole feature is presented in Fig- 
ure [3] The hole is classified by Boomsma (2007) as Type 1 under the scheme described 



by Brinks & Bajaja (1986), meaning that it is characterized by a gap in the PV diagram 
but without the appearance of an expanding bubble. This particular hole is most easily 
recognized in Figure [3] by inspection of the contours and noting the density enhancement at 
the edges of the hole. A two-component fit to the average velocity profile over the width of 
the hole is also shown. The combination of two Gaussian components produces an excellent 
match to the data. The velocity difference between the components is 11.7kms~^; correct- 
ing for the inclination angle of 38°, the velocity along the axis perpendicular to the disk 
is 14.8 km s~^. The velocity dispersions of the two components increases in the redshifted 
component to 14.1 km s~^ compared to 6.5kms~^ in lower- velocity component, which we 
associate with the ISM in the main disk. Assuming that the redshifted gas is located on the 
front side along with the magnetic feature, the second H I component represents gas which is 
falling back onto the disk at an average velocity of 15kms~^. Inspection of the PV diagram 
shows that the infall is more pronounced on the SE (radially outward) side of the hole. Such 
wings in the velocity profiles are ubiquitous in NGC 6946 and are certainly not unique to 
this region, but this particular redshifted feature fits naturally into the picture described 
below. Note that a blueshifted component, if present, would be confused by Hi emission 
from the Milky Way as indicated in Figure |3} 

The magnetic field strength in the vicinity of the bubble can be estimated by taking 
advantage of the fact that the observed RM gradient (ARM) is symmetric about the average 
RM in the region. Half of the RM gradient can thereby be attributed to the field strength, 
which is to say that 

^ ARM = 0.81 / n,B -dL (1] 



2 

Taking the inclination i of the galaxy into account, we have 

^ ARM = 0.81ne |5| L sinasinz, (2) 

where L is the path length through the volume of interest and a is the angle between the 
magnetic field lines and the disk plane. Taking the inclination to be 38°, and estimating 
a = 45°, L = 300 pc (half of the distance across the RM feature), and a typical thermal 



electron density Ue = 0.05 cm (cf. Ferriere 2001), the resulting regular magnetic field 
strength is determined to be 
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Major axis offset (kpc) Velocity (km/s) 

Fig. 3. — Kinematics of the H I hole. Left: PV diagram along the slice shown in the top left 
panel of Figure |2} The orientation along the slice (SE to NW) is indicated. Contours start 
at 0.36mJybeam~^ (l.Sa) and increase by powers of two. Dotted contours are drawn for 
negative values. The locations of the edge of the H I hole are indicated with vertical dashed 
blue lines. Note the high-velocity wing at the location of the hole. The red hatched region 
(at negative velocities) shows the area of the diagram which is contaminated by MW H I 
emission. Right: the average velocity profile within the hole region of the PV diagram, along 
with a two-component Gaussian fit to the profile. The red and blue components together 
give the combined fit, which is shown in green. The residuals of the fit are shown at the 
bottom. Again, the red-hatched region indicates the values that are contaminated by MW 
Hi (that region was not included in the fit). 
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which compares well with the magnetic field strength previously measured in the range 



8 — 10 /iG in the magnetic arms (Beck 2007), as is the feature studied here. 



Taken together, these lines of evidence provide new observational support to a picture 



consistent with the chimney model (Norman & Ikeuchi 1989) of disk-halo interaction, in 



which star formation activity in the disk drives convective hot gas fiows upward into the 
halo, carrying magnetic fields along with the hot gas motion, and forming bubbles in the 
disk that are observable as H I holes. The hot gas cools in the halo, condensing into H I clouds 
which then return to the disk on time scales of a few tens of Myr and with a net outward 



radial movement (Collins et al. 2002 Fraternali & Binney 2006). State-of-the-art numerical 



models (e.g., Korpi et al. 1999 de Avillez fc Breitschwerdt|2005 ) depict the detailed evolution 
of the ISM during this process. The observational results presented here make the first clear 
identification of magnetic involvement in disk-halo fiows. This has significant consequences 
for our understanding of the galactic dynamo process, because large-scale magnetic fields 
are carried upward by the hot gas, but will not return downward with the cool clouds 
traced by the redshifted H I ( Brandenburg et al.|1995 ). Moreover, the smaller-scale turbulent 
component of the magnetic field will also be carried out of the disk, which has been proposed 



as a mechanism to prevent quenching of the dynamo process (Shukurov et al. 2006). As a 



related point, we remark that a close connection in NGC 6946 between the position angle of 
the elliptically shaped H I holes and the local spiral pitch angle has previously been noticed 



(Boomsma et al. 2008), and may actually refiect the fact that superbubbles preferentially 



expand along the direction of the ordered field (Tomisaka 1998 Stil et al. 2009) 



It is remarkable that of the many Hi holes in NGC 6946 (121 are listed by Boomsma 



et al. 2008), only one shows a clear RM signature. In estimating how unique this feature 



may be, there are several factors to take into consideration. 



The lack of polarized synchrotron radiation on the receding (southwest) side of the 
galaxy, clearly seen in Figure [l| is now understood to be a geometric effect (Braun 
et al. 2010), so no RM information is available for the subset of holes in that quadrant 



of the disk. This effect is relevant for about 30 (25%) of the catalogued holes. 

Many holes are located where more intense star formation is taking place, which in- 
creases the level of turbulence and thus the local degree of depolarization. This explains 
why there is an overall anticorrelation between the locations of H I holes and regions 
with diffuse polarized synchrotron radiation. It is likely that the particular feature 
studied here is located some distance above (closer to the observer than) the turbulent 
midplane of NGC 6946, making it visible in diffuse polarization. Only 16 (13%) of the 
catalogued H I holes fall fully within the regions which do not suffer from depolariza- 
tion effects or the geometric effect mentioned above (see Figure [2]). Estimating that 
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the diffuse polarized emission originates at least one exponential scale height above the 
midplane on the frontside, only about 25% of those 16 holes could have been traced 
by diffuse polarized emission. 

The age of the hole must be within a certain range: old enough that the field con- 
figuration has had time to gain a significant vertical offset, but young enough that 
vertical shear has not yet destroyed the observational signature. For an outflow speed 
of lOOkms"^, the feature can grow to 300 pc height after only about 3Myr. The ver- 



tical shear is not known for NGC 6946 (Boomsma et al. 2008), but it is suspected to 



be at least as high as the measured value in NGC 891, ~ 15kms ^ kpc ^ (Heald et al. 



2006 Oosterloo et al. 2007). For this, the characteristic shear time would be around 



60Myr. Other galaxies have higher shear values (e.g., Heald et al.||2007 ), so the length 
of the observational window is short indeed but matches well with the age estimated 
for this particular hole from its H I properties. 

Finally, we note that the azimuthal location of the H I hole within the disk of NGC 6946 
is only 12° from the minor axis. Since the pitch angle of the spiral pattern in the 



magnetic field vectors is about 18° (cf. the optical spiral pitch angle given by Kennicutt 



(1981), 28°), this means that the location in the disk leads to the ideal geometrical 



situation for identifying an RM gradient across the feature. At other locations in the 
disk such a feature would have a different and far less obvious observational signature. 



Considering all of these effects, we estimate that the feature described here could have been 
detected in only about 4 of the 121 catalogued H I holes. It therefore seems reasonable to 
expect that the vertical magnetic field transport detected in this one H I hole is actually 
much more ubiquitous across the disk of NGC 6946 and other similar galaxies. Sensitive 
observations at higher radio frequencies, where depolarization effects are less severe, can be 
utilized to test this prediction. Observations in the 2 — 4 GHz range, for example, would 
provide a reasonable compromise between recovery of depolarized diffuse emission on the 
one hand, and acceptable RM precision on the other hand. 
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